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Y. K. Ho and Ronald J. W. Henry 
Department of Physics and Astronomy 
Louisiana State University 
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Abstract 

Oscillator strengths between various doublet states of Oil 
ions are calculated in which extensive multi-conf iguration wave 
functions are used. The lower levels for the transitions are of 
the 2p* Z D° and 2p 3 2 P° states, and the upper levels are 2p“ , 3s, 
and 3d states. The results, which are estimated to have errors of 
less than 10% for individual transitions, agree quite well with 
the beam foil experiments, as well as with the calculations by 
use of the non-*closed shell many electron theory (NCMET). The 

agreement with the rocket measurements is also good except for 

• e 

the 538/581 A pair, in which the 538 A line is believed to be 
blend with the other stronger quartet line. However, a compari- 


son with the recent 
discrepances between 
do exist for certain 


branching ratio measurement 
the present calculation and 
transistions. 


indicates that 
the experiment 


Cl) Introduction 


We report a calculation for the oscillator strengths of the 
Oil ions between various doublet states by the use of extensive 
configuration interaction wave functions. Studies of various 
oxygen ions have practical, experimental, and theoretical intei — 
ests. From the practical side, accurate atomic data are essen- 
tial for the interpretations of the physical conditions of astro- 
physical plasmas. For example, the recent Voyager flyby of the 
Jupiter has revealed bright extreme-ultraviolet (EUV) lines of 
oxygen and sulfur ions in the vicinity of the satellite Io 1 . In 
addition, some EUV transistions between various doublet states of 
Oil ions have been identified in the day airgrow spectrum. The 
spectrum was observed in a rocket measurement 2 , 3 . From the 
experimental side the use of beam foil techniques to determine 
the lifetimes of variois ionic states has been quite popular 41 , 5 . 
Other recent experimental investigations on the Oil ions include 
the branching ratio measurements*, 7 . On the theoretical side, 
accurate calculations of the oscillator strengths for Oil ions 
are of recent interest. Since OH ions are singly ionized spec- 
ies, strong configuration interactions exist in the calculations 
of oscillator strengths. Furthermore, since Oil are also open 
shell ions, various correlation effects must be included for 
accurate calculations. Previous theoretical investigations 
include the calculations by use of the non-closed shell many 
electron theory (NCMET) 3 and the first order theory of oscillator 
strengths (F0T0S) 9 , which is itself a subset of the full NCMET 
theory. The oscillator strengths for these open shell and singly 
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ionized ions are sensitive to the wave functions used, as well as 
to the procedures by which the parameters in the wave functions 
are determined (optimized). Such sensitivities have been docu- 
mented in the literature 1 0 , 1 1 , 1 2 . The results for some of the 
transitions obtained by these two calculations differ from each 
other considerably. It is felt, because of these experimental 
and theoretical interests, that an independent extensive calcula- 
tion for the oscillator strengths of Oil ions is worthwhile. 

This work presents such a calculation. The wave functions 
used are of extensive conf igurat ion-interact ion type. The pro- 
gram CIV3 of Hibbert 13 is used in the present calculation. In 
section (II) we will describe the wave functions as well as the 
procedures to determine the parameters for the wave functions. 
In section (III) we will present our results. Comparisons with 
other calculations and with experiments will also be given. 

(II) Theory and Calculations 

In the program CIV3 of Hibbert 13 , the wave functions are 
expressed as linear combinations of configuration-interaction- 
type wave functions. 


* = 2 a, i, , 


( 1 ) 


with constructed from one-electron orbitals 
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-1 


( 2 ) 


u = r P (r) Y (6,^) x(m ) . 
nlm nl lm s 


In Eq.(2), u is a product of a radial function, a spherical hai — 
monic and a spin function. The radial functions are linear com- 
binations of Slater orbitals 


k I j -«jr 

P (r) = 2 C r e (3) 

nl j=1 jnl 


Uave functions are obtained by solving the eigenvalue prob- 
lem 


<4'|H|4’> = $> E (4: 

i j ij j 

We will discuss later in this article how various parameters for 
the orbital wave functions are determined. Once the wave func- 
tions for both the upper and lower states are obtained, the 
absorption oscillator strengths can be calculated. The oscilla- 
tor strengths (expressed in atomic units) in the length and 
velocity forms are; 


AE 2 

f = 2/3 | <4- |Zk r | 4 1 >| (length) (5) 

L g i k f 

i 


( aE )- 1 _ 2 

f = 2/3 | <4- |2 k p >| (velocity) (6) 

V g i k f 

i 


where r and p are the coordinate and momentum operators 
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respectively, and aE the energy difference between the initial 

(+ ) and and final ( + ) states. The summations in equations (5) 
i f 

and (6) are summed over the number of electrons. A necessary 
(but not sufficient) condition for an accurate calculation for 
the oscillator strengths is the consistence between the length 
and velocity values. 

The procedure to choose the configurations to represent vai — 
ious ionic states is the following; Since Oil is an open shell 
ion, we include the three major correlation effects in the pres- 
ent calculation. These effects are; 

(1) internal effects that include the configurations for 
which the electrons occupy the Hartree-Fock (H-F) sea, (For exam- 
ple, for the 2p 3 2 P° state the H-F sea includes the Is, 2s and 2 p 
orbitals. Therefore, the 2p 5 2 P° configuration represents the 
internal correlation effect.) 

(2) semi-internal effects include the conf igurat ions in 
which only one electron is outside the H-F sea, and 

(3) all-external effects include the configurations in which 
two or more electrons are outside the H-F sea. 

In this work, however, we limit the all-external effects to 
the conf igurat ions in which no more than two electrons are out- 
side the H-F sea. Such an limitation is called the quasi-extei — 
nal effect by Sinanoglu" . 

For illustrative purposes, we show in table 1 the configura- 
tions that are used to represent the 2S2P 4 * 2 D state together with 
their classif ications. For the semi-internal correlation effects 
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we include contributions from both the N=3 shell and the N=4 
shell. However, for the quasi-external effects only the contri- 
butions from the N=3 shell (except for 4f) are included. Also 
implicitly, we do have some contributions from higher shells. 
Since some of the orbitals do not necessarily have spectroscopic 
meanings, the '31' orbital therefore represents the average con- 
tributions from the higher shells. 

The parameters for different orbitals are optimized individu- 
ally according to different transitions. For example, for the 
transition between the 2 p 3 2 D° - 2 D states, we use the Is, 2s and 
2p orbitals that are optimized on the 2p* 2 D C state in the single 
configuration H-F approximation. The parameters for these orbi- 
tals have been published by Clementi and Roetti 1 *. The 3s orbi- 
tal is then optimized on the 3s 2 D excited state, and the 3p and 
3d orbitals, etc., are optimized on the 2p* 2 D excited state. 
Therefore, in this case, only the 3s orbital has spectroscopic 
meaning, while others are simply correlation orbitals. 

In order to avoid the variational collapse*, care must be 
taken for a state that is not the lowest state in a given symme- 
try (the totai angular momentum, total spin, and parity). For 
example, in the optimization of the 3s orbital for the 2p 2 3s 2 D 
state (which has the second lowest energy in its symmetry) , the 
configuration of the 2s2p u is also included. We then optimize 
the 3s on the second lowest energy eigenstate. The Hyl leraas-Un- 
dheim-MacDonald theorem would guarantee that the calculated 
energy be an upper bound to the true energy for the 3s 2 D state. 
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For the calculations of the jd states, we optimized various 
orbitals by the following procedure. Since 3d also contributes 
to the ground doublet states, we optimize the 3d orbitals on the 
average energy of the ground state and the 3d excited state. The 
4d orbital was then optimized on the 3d excited state. As a 
result, the 3d excited state is represented by the combination cf 
3d and 4d. Here again the individual nd orbital does not have 
spectroscopic meanings. Furthermore, in the optimization proce- 
dure we also include the 3s and the 2p* states in the expansion. 
Since the 3d state has the third lowest energy for a given set of 
(tt.L.S), we have to include the two lower states to avoid the 
possibilities of variational collapse. 

Most of the transitions reported in this work are done by the 
procedures described above. In other words we first construct a 
good representation for the lower state. The upper state is rep- 
resented by the extensive conf igurat ion-interact ion type wave 
functions that include the internal, semi-internal and quasi-ex- 
ternal effects. In such a procedure we have found the theoreti- 
cal energy difference agrees quite well with the experimental 
observed wave length for most of the transitions reported in this 
work. However, for the transition between the 2 P° - 2 S states we 
have found that the theoretical energy differences between the 
upper and lower states are not consistent with the experimental 
values. An explanation is suggested as the following; Since the 
2 S states are located at higher positions, the orbitals which are 
good for the lower 2p* 2 P° state may not be good for the upper 
states. Therefore we treat the 2p 2 3s 2 S state as a state with 
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an 3s electron attaching to the 2p* ‘S core of the OIII ion. As 
a result, the orbital parameters for the Is, 2s and 2p are those 
of the OIII ‘S state. The 3s orbital is then optimized on the 3s 
*S excited state. By doing this we have found that the energy 

differences between the lower and upper states are consistent 
with the experimental wave lengths. 

(Ill) Results and Discussions 

We present our results in table 2 as well as those by Sina- 

roglu and coworkers who used the program "ATOM”. Such a program 

is based on the non-closed shell many electron theory (NCMET)*. 

The general agreements between the present calculations and the 

f values of the NCMET results are quite good, with tne present 
L 

f and f values being more consistent than those in Ref. (Cj. 
L V 

It should ^e mentioned that these two sets of results are 
obtained by two different programs, although the methods of 
selecting various conf igurations are very similar. Also, by 
judging the differences between the oscillator strengths in the 
length and velocity forms, we estimate the errors for the present 
calculations are within 10?! for the individual transitions. 

The comparisions with the F0T0S* calculations are good except 

© e 

for the 673 A and 538 A transitions. In these cases substantial 
disagreements have been found. This may reflect the sensitivity 
of the optimization procedures (for the various orbital parame- 
ters) used in Ref. (9). It should also be mentioned that some of 
the transitions reported in this work have been calculated by the 
present authors using less extensive wave functions 18 . Also in 
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R«f . (15) the parameters used to represent variojs state wave 
functions were optimized for the quartet states. The present 
results hence represent improvements over the previous calcula- 
tions in Ref. ( 15) . 

Comparisions of the present results with different experi- 
ments are shown in tables II and III. It is seen that the pres- 
ent results agree quite well with those by the beam foil experi- 
ments 1 (table II). The comparision with the branching ratio 
experiment (table III) is less satisfactory for certain tran- 
sitions. Since the errors for the present calculations are esti- 
mated to be within 105 for the individual transitions, the error 
for a given ratio pair is hence of 205. It is seen that some of 
the present results differ from the branching ratio experiment 7 , 
even when the combined experimental and theoretical errors are 
included. For example, the discrepancy for the 482/ 515 A pair 
shows a factor of two difference, and the 555/601 A pair shows a 
difference of about 605. Of course we also have good agreement 

9 

for the 538/581 A pair. Generally the present ratios (with the 
longer wavelength as the denominator) are smaller than the 
branching ratio experiments. However, the present results seem 
to agree better with the EUV airglow rocket measurements, except 

o s 

for the 538/581 A ratio. The observed 538 A line (2p* 2 P - 2p* 
*D°) is believed to be blend with the 538 A (3s *P - 2p* ‘S 0 ' 
line where the latter has a larger oscillator strength. 

This work presents an extensive calculation for the oscilla- 
tor strengths for Oil ions between various doublet states. The 
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discrepancies between the theoretical calculations and the 
branching ratio experiments indicate that further theoretical and 
experimental investigations on such systems are worthwhile. From 
the theoretical point of view, improvements of the wave functions 
may be made if one also includes (1) the semi-internal correla- 
tion effects coming from the higher shells 'higher than N=4), ( 2 t 
the external effects for which more than two electrons are out- 
side the H-F sea, and (3) the quasi-external contributions coming 
from higher shells (higher than N*3). Of course, such improve- 
ments may have small effects on the oscillator strengths. It is 
hoped, therefore, that the present calculations would stimulate 
further experimental investigations, which in turn would provide 
a more stringent criterion for the test of various theoretical 
calculations . 
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TABLE I - Configurations used to represent the 2s2p D state 

Internal 

2s2p 4 


Semi -internal 


2s 2 2p 2 3s 

2s 2 2p 2 4s 

2p 4 3s 

2p 4 4s 

2s 2 2p 2 3d 

2s 2 2p 2 4d 

2p 4 3d 

2p 4 4d 

2s 2p 3 3p 

2s2p 3 4f 

Quasi-external 

2s2p 2 3s 2 

2s2p 2 3s3d 

2s2p 2 3p 2 

2s2p 2 3p4f 

2s2p 2 3d 2 

2p 3 3s4f 

2s2p 2 4f 2 

2s 2 2p3s3p 

2p 3 3s3p 

2s 2 2p3p3d 

2p 3 3p3d 

2s 2 2p3d4f 

2p 3 3d4f 

2s 2 2n3s4f 
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TABLE III Comparisons with experimental measurements 
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